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Abstract 

We study the production of 7 rays via LSP annihilations in the core of the Galaxy as a possible 
experimental signature of the constrained minimal super symmetric extension of the Standard Model 
(CMSSM), in which super symmetry-breaking parameters are assumed to be universal at the GUT 
scale, assuming also that the LSP is the lightest neutralino x- The part of the CMSSM parameter 
space that is compatible with the measured astrophysical density of cold dark matter is known 
to include a fi — x coannihilation strip, a focus-point strip where x enhanced Higgsino 

component, and a funnel at large tan/3 where the annihilation rate is enhanced by the poles of 
nearby heavy MSSM Higgs bosons, A/H. We calculate the total annihilation rates, the fractions 
of annihilations into different Standard Model final states and the resulting fluxes of 7 rays for 
CMSSM scenarios along these strips. We observe that typical annihilation rates are much smaller 
in the coannihilation strip for tan/3 = 10 than along the focus-point strip or for tan /3 = 55, and that 
the annihilation branching ratios differ greatly between the different dark matter strips. Whereas 
the current Fermi-LAT data are not sensitive to any of the CMSSM scenarios studied, and the 
calculated 7-ray fluxes are probably unobservably low along the coannihilation strip for tan /3 = 10, 
we find that substantial portions of the focus-point strips and rapid-annihilation funnel regions 
could be pressured by several more years of Fermi-LAT data, if understanding of the astrophysical 
background and/or systematic uncertainties can be improved in parallel. 
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1 Introduction 



Relatively soon after the realization that the lightest supersymmetric particle (LSP), in par- 
ticular the lightest neutralino Xi could naturally provide the astrophysical dark matter in 
models in which i?-parity is conserved [1], it was suggested that LSP dark matter annihila- 
tions might be detectable via features in the 7-ray spectrum [2]. This is now a hot theoreti- 
cal |3]-[TTj and experimental topic, with several experiments [T214T9] studying the cosmic-ray 
7 spectra from a variety of astrophysical sources such as the galactic centre and bulge as well 
as dwarf galaxies. There have been claims of deviations from calculations of conventional 7- 
ray backgrounds, and corresponding claims of evidence for new physics invoking tailor-made 
supersymmetric scenarios [20]. In parallel, particularly in the context of searches with ac- 
celerator experiments, there have been many studies of simplified supersymmetric scenarios. 
Foremost among such scenarios is the constrained minimal supersymmetric extension of the 
Standard Model (CMSSM), in which the soft supersymmetry-breaking parameters mg, 7^11/2 
and Aq are assumed to be universal at the supersymmetric GUT scale [2TH26] . 

Accelerator experiments and astrophysical searches have complementary roles to play 
in elucidating the nature of any dark matter particles. Accelerator experiments cannot 
determine whether any candidate dark matter particle is in fact stable, rather than merely 
living long enough to escape from the apparatus, while astroparticle experiments are limited 
in their capabilities to disentangle the dynamics of dark matter models and thereby, for 
example, verify that they yield the appropriate cosmological dark matter density. Connecting 
the accelerator and astroparticle experiments requires interpreting them within a common 
model framework, and in this paper we choose to compare their sensitivities within the 
framework of the CMSSM, assuming that i?-parity is conserved and that the LSP is the 
lightest neutralino x- 

The parameters of the CMSSM are the universal scalar mass mo, gaugino mass mi/2 and 
trilinear supersymmetry-breaking parameter Aq, the ratio of Higgs vevs tan/3 and the sign 
of the Higgs mixing parameter /x. Here we focus on /i > 0, motivated by — 2 [271[28] and 
(to a lesser extent) by 6 — j- S7 [29]. We study the cases tan/3 = 10,55, which bracket the 
phenomenologically-plausible range. Currently, there is little experimental sensitivity to Aq, 
and we limit our discussion here to the case = 0. 

As is well known, the regions of the CMSSM parameter space in which the dark matter 
density falls within the narrow range permitted by WMAP and other cosmological observa- 
tions [30] may be represented as relatively narrow strips in (mi/2, frio) planes for fixed values 
of tan /3 [25|[3T] , which are illustrated in Fig. [TJ As seen in the left panel, if tan /3 = 10 there 
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is one WMAP strip with mo ^ fni/2 where the rehc density is reduced into the allowed 
range by coannihilations of the relic neutralinos x with sleptons i [32] , and another strip at 
relatively large tuq in the so-called 'focus-point' region [33] where the relic density is brought 
into the allowed range by enhanced annihilation due to a relatively large Higgsino component 
in the composition of the lightest neutralino x- If t^ii P = 55, as seen in the right panel of 
Fig. [H the coannihilation strip segues into a funnel region [2T1|23] where the relic density is 
brought into the allowed range by rapid annihilations through direct-channel heavy Higgs 
resonances H/A. Also visible in Fig. [1] are the constraints imposed by 6 — s'-f, the LEP lower 
limits on chargino and Higgs masses [3l], and the region favoured hj — 2 if the apparent 
experimental discrepancy with the Standard Model calculation is ascribed to supersymme- 
try. Fig. [1] also shows the constraints implied by the absence of any supersymmetric signal 
in the 2010 LHC data [35]. 

In this paper we study the likely sensitivity of searches for 7 rays from the galactic centre 
to LSP annihilations in the CMSSM, taking as examples the (mi/2, fno) planes shown in Fig.[T] 
and focusing on the dark blue WMAP strips, in particular. The first step is to understand 
relevant features of CMSSM dark matter annihilation processes, which we study in Section 2. 
As we discuss there, the x~X annihilation cross section is much smaller in the coannihilation 
region of the CMSSM for tan/3 = 10 than it is in the focus-point region, or for tan/3 = 55. 
This is easily understood, because in the tan /3 = 10 coannihilation strip the relic density 
is brought down into the WMAP range by the sum over all x — £ coannihilation processes, 
relative to which x~X annihilation is numerically small. Moreover, only S-wave annihilation 
is important in the Universe today, whereas P-wave annihilation also played a role in the early 
Universe, particularly for the parameter space associated with the co-annihilation strip. The 
smallness of the S-wave x~X annihilation cross section in the tan /3 = 10 coannihilation strip 
implies that all astroparticle searches for dark matter annihilation products, whether they 
be photons, neutrinos, positrons, antiprotons or antideuterons, will be relatively insensitive 
in this region of the CMSSM parameter space. 

A second issue we study in Section 2 is that of the branching fractions for annihilations 
into different final Standard Model states. Indicative explorations of annihilation signatures 
may assume particular final states as illustrations, but in the CMSSM these branching frac- 
tions are fixed at each point in parameter space. We find that the dominant annihilation 
channels are r+r^ in the coannihilation region and W^W~ in the focus-point region for 
tan /3 = 10, and hh in both regions for tan /3 = 55. The 77 final state, which would be 
particularly promising for detection via a 7 line with = rriy, unfortunately has a very 
small branching fraction in all cases. The next step is to model the hadronic, leptonic and 
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Figure 1: The {1^1/2,1^0) planes in the CMSSM for ta.n f3 = 10 (left) and tan P = 55 (right), 
assuming /i > and Aq = 0, showing the 95% CL limits imposed by ATLAS and CMS data 
(purple and yellow lines, respectively). The regions where the relic LSP density falls within 
the range allowed by WMAP and other cosmological observations appear as strips shaded 
dark blue. The constraints due to the absences of charginos and the Higgs boson at LEP 
are also shown, as black dashed and red dot-dashed lines, respectively. Regions excluded by 
the requirements of electroweak (EW) symmetry breaking and a neutral LSP are shaded dark 
pink and brown, respectively. The green region is excluded by b ^ s'~f, and the pink region 
is favoured by the supersymmetric interpretation of the discrepancy between the Standard 
Model calculation and the experimental measurement of — 2 within ±1 and ±2 standard 
deviations (dashed and solid lines, respectively). 



7 components of the final states, for which we use PYTHIA, as introduced at the end of 
Section 2. 

In Section 3, we turn to the treatment of the astrophysical background and astrophys- 
ical aspects of indirect detection. In particular, we discuss possible dark matter profiles 
towards the center of the galaxy, including Navarro-Frenk- White (NEW) [36] and Einasto 
models [371[38], as well as a simple isothermal model [39]. We also compare our predicted 
signal with a model background and the current data and possible sensitivity of the Eermi- 
LAT detector [13] in Section 3, illustrating our discussion with some specific benchmark 
CMSSM scenarios. Then, in Section 4, we present estimates of the perspectives Eermi-LAT 
may have for distinguishing a possible annihilation 7-ray signature from the background 
along the WMAP strips shown in Eig. [H As might be expected on the basis of the magni- 
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tudes of the annihilation cross sections in the different regions of CMSSM parameter space, 
for tan /3 = 10 there are better prospects for disentanghng a signal from the background in 
the focus-point region, whereas for tan /3 = 55 there are prospects in both the funnel and 
focus-point regions. However, putting pressure on the CMSSM scenarios studied here with 
future years of Fermi-LAT data would require considerable improvements in the modelling 
of the background and/or reduction in the systematic errors. Section 5 summarizes our con- 
clusions about searches for 7 rays from the galactic centre, and makes some remarks about 
other search strategies for CMSSM dark matter annihilations. 

2 CMSSM Dark Matter Annihilation Processes 

The relic density of neutralinos is determined by annihilations when they are slightly non- 
relativistic. Typically, annihilation freeze-out occurs when x = T/m^ ~ "^rci/S ~ 1/23, 
where Vj-e\ is the relative velocity of annihilating particles and T is the temperature of the 
Universe. For small x < 1, the annihilation cross section can be expanded in a series as 
(crWrci) = a + bx + . . . [HllQ]. Because neutralinos are Majorana particles, we generally have 
a oc m'l/mt,^^,, where rrif is the mass of a final-state fermion, and b oc m^Jmt,^^,. Hence, often 
a <^b and the relic density is largely determined by 6, with higher-order coefficients in the 
expansion relatively unimportant. Any signal of dark matter annihilations will also be scaled 
by the annihilation cross section. However, since dark matter particles annihilating in the 
galactic halo today are very non-relativistic with fj-ci ^ 1, the annihilations are essentially 
pure S-wave, and (afrci) ~ o- This is a valid argument even close to the inner region of the 
galactic halo, where the N-body simulations indicate that the DM halo particles are quite 
non-relativistic, with fj-ci ~ 0(10"'^ — 10""^) [36] . 

We note in passing that in our calculation the Sommerfeld enhancement effect [S] is 
not important. One condition for having such a sizeable enhancement is a high degree 
of mass degeneracy between the LSP and another sparticle, such as the chargino or stau. 
However, along the focus-point WMAP strip it is not possible to achieve the required amount 
of degeneracy between the LSP and lightest chargino and, on the other hand, along the 
stau coannihilation strip the effect of the Sommerfeld enhancement is known not to be 
important [12]. 

With this in mind, we display in Fig. [2] the S-wave x~X annihilation cross section along 
the WMAP-compatible strips shown in Fig. [H as a function of mi/2 in each case. The 
coannihilation/funnel strips are represented by solid lines, black for tan/3 = 10 and red for 
tan /3 = 55, and the focus-point strips are represented by a blue dotted line for tan /3 = 10 
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Figure 2: The x^X annihilation cross section along the WMAP strips in the coannihilation, 
focus-point and funnel regions for tan /3 = 10,55, = and fi > 0, as functions of mi/2- 
We see that the annihilation cross section along the tan /3 = 10 coannihilation strip is much 
smaller than along the other strips, and decreases rapidly as mi/2 increases. 



and a mauve dashed line for tan/3 = 55. In each case, we have fixed Aq = and for each 
value of mi/2 we have adjusted the value of mo to obtain the WMAP value of the relic 
density, n^h"^ = 0.1109 ± 0.0056 [30]. 

We notice immediately that the total x ~ X annihilation cross section is much smaller 
along the coannihilation strip for tan/3 = 10 than in the other cases, falling rapidly as mi/2 
increases. This refiects the importance of x ~ ^ coannihilations in reducing the cosmological 
relic density into the WMAP-compatible range, a role that becomes increasingly important 
at larger mi/2 where the x ^"^^ ^'^'s increasingly degenerate. We recall that, as seen in 
Fig. [H the LHC upper limits on sparticle production enforce mi/2 > 350 GeV along this 
strip, where av < 10~^^ cm^/s, whereas — 2 favours the portions of the coannihilation 
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strips with mi/2 < 400 GeV for tan/3 = 10 and mi/2 < 800 GeV for tan/3 = 55. 



ong the 
This 



On the other hand, the total S-wave annihilation cross section is much larger aL 
other WMAP strips shown in Fig. [2], with typically crt" ~ (1 — 2) x 10^^^ cm^/s 
difference immediately suggests that detecting dark matter annihilations will be easier along 
the funnel for tan/3 = 55 as well as the focus-point strip for tan/3 = 10, or for tan/3 = 55. 
We note, however, that these regions are disfavoured by fits to — 2 and other low-energy 
precision observables [13]. 

In the case of the coannihilation/funnel strip for tan/3 = 55, we note the appearance of 
a second red line for mi/2 > 1200 GeV. This reflects the appearance in Fig. [2] of a second 
branch of the WMAP strip on the other side of the rapid-annihilation H /A funnel seen 
in the right panel of Fig. [H The annihilation cross section takes similar values along the 
focus-point strips for both the tan/3 = 10 and 55 cases (except at small mi/2), and the LHC 
limits have no impact along either of these strips. As seen in Fig. ^ — 2 favours only the 
portion of the tan /3 = 55 focus-point strip with mi/2 < 200 GeV, which is disfavoured by 
other constraints [^, and disfavours all the tan /3 = 10 focus-point strip. 

The detectability oi x ~ X annihilation depends also on the branching fractions for an- 
nihilations into specific Standard Model final states and the 7 spectra they produce. Fig. [3] 
displays the branching fractions for the most important final states, and we see that they 
are quite different along the various WMAP strips studied. In the case of the coannihilation 
strip for tan /3 = 10 (upper left panel), we see that r+r^ final states dominate at low mi/2, 
followed by hh final states, with W^W~ and it final states gaining in importance at larger 
mi/2, where the total annihilation cross section is, however, much reduced as seen in Fig. [2l 

On the other hand, the roles of the hh and r"*"r~ final states are largely reversed along the 
coannihilation/funnel strip for tan/3 = 55 (upper right panel), including the second branch 
when mi/2 ~ 1500 GeV. In the case of the focus-point strip for tan /3 = 10 (lower left panel), 
we see that W^W~ final states dominate, except for hh final states at small mi/2 ~ 200 GeV 
and the appearance of it final states when mi/2 > 400 GeV. Finally, in the case of the 
focus-point strip for tan/3 = 55 (lower right panel), we see that hh final states dominate over 
r+r~ final states everywhere. 

The prompt component of the photon flux, i.e., the photons that are produced directly 
by neutralino pair annihilation in the halo, consists of two components: the monochromatic 
and the continuum. The monochromatic component due to the 77 branching fraction is very 



^There is an exception at small mi/2 on the tan/3 = 10 focus-point strip, where the annihilation xx ^ 
W^W~ is kinematically inaccessible. 

^Wc note that, for tan/3 = 55, the LHC limits have less impact on mi/2 than the 6—5-57 constraint, 
which imposes toi/2 > 400 GeV. 
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Figure 3: The branching fractions for x ~ X annihilations into pairs of Standard Model 
particles along the WMAP strips for tan/3 = 10 (left panels) and tan/3 = 55 (right panels), 
in the coannihilation and funnel regions (upper panels) and in the focus-point region (lower 
panels). 



small in all cases. For tan /3 = 10, the branching fraction is less then 10^'^ for < 400 
GeV. It rises slowly to roughly 0.015 as the endpoint of the coannihilation strip is reached, 
but here the total annihilation rate is very small, as seen in Fig. El In the focus-point region 
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for tan/3 = 10, the branching ratio is always neghgible and remains below 10"''. In the 
case of the coannihilation/funnel strip for tan/3 = 55 it remains below a few xlO~^, and 
is even lower along the focus-point strip for tan /3 = 55. Further details, especially for the 
monochromatic part from the processes XX ~^ '^l/^li be given in the next section. 

The continuum component is due to photons produced by the multibody neutralino pair 
annihilation channels, and is usually modelled using event generators such as PYTHIA [S] or 
HERWIG [l5]. Here we use PYTHIA to calculate the continuum component of the differential 
photon flux dN^^^I dE^. The decay channels we take into account are: qq for all quark flavors, 
and the r+r", W~^W~, ZZ, gg, Zh and ZH flnal states B In particular, we follow [16] and 
references therein, parametrizing the flts of photon spectra produced by PYTHIA, in order to 
obtain the differential flux as an analytic function of the parameter x = E^/m^, for each LSP 
mass m^. The functions dN^^^^/dx for each flnal state are related to the energy spectrum 
used below dN^^^^/ dE^ = {1 / m^)dNl^^^/ dx . Our results for the continuum component agree 
very well with those of [46ll47j. In the modeled fluxes we have included effects that are 
related to flnal state radiation [4]. Recently, it has been shown that flnal state strong and 
EW corrections from the radiative emission of gluons, W or Z |^ can enhance otherwise 
P-wave suppressed channels, like the neutralino pair annihilation to electrons and neutrinos. 
In principle, corrections to the neutralino annihilation to electrons can affect the photon 
flux, due to the Inverse Compton Scattering (ICS) effect, but as we will see in the following 
section, this part of the signal is subdominant for energetic photons. 

3 Gamma fluxes and the Fermi-LAT data analysis 

In order to estimate the sensitivity of a detector like Fermi-LAT [13] to new sources of photon 
fluxes arriving from the galactic center (GC), one must calculate both the signal and the 
background contributions to the photon flux, which we do in the following subsections. 

The GC is probably one of the most complicated regions of the Milky Way, from both the 
observational and the theoretical points of view. Observationally, the emitted electromag- 
netic spectrum exhibits various interesting features, like the recently-observed giant 7-ray 
bubbles, the so-called "Fermi bubbles" [49], that can be related to accretion activity. From 
the theoretical point of view, the supposed existence of a massive black hole in the GC might 
change drastically the size of the 7-ray flux [7]. Moreover, other phenomena such as turbu- 
lent galactic winds, large magnetic flelds etc, can affect the theoretical predictions [50]. In 

■^In the regions of the CMSSM parameter space we study here, the ZH final state, where H is the heavy 
CP-cven Higgs boson, yields a very small contribution. 
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this work we adopt a conservative approach. In particular, we do not consider the possible 
effects due to the presence of a massive black hole, or any other source that might enhance 
the 7-ray flux arising from annihilations of neutralinos in the dark matter halo. 



3.1 The supersymmetric signal 

In the context of supersymmetric models where the neutralino plays the role of the dark 
matter particle, there are two main components in the signal: 

• Prompt photons produced either directly through one-loop processes such as 2% — )■ 
27/Z7 (monochromatic components) or indirectly as final state radiation or through 
the hadronization, fragmentation and decays of the tree-level neutralino pair annihila- 
tion products discussed in the previous section (the continuum component); 

• Photons that originate from energetic electrons and positrons produced in neutralino 
pair annihilation via the Inverse Compton Scattering (ICS) process on background 
photons. 

There is also a synchrotron component produced by the energetic electrons and positrons 
accelerated in the galactic magnetic field, but it contributes at energies much smaller than 
the prompt and ICS components [51], below the effective Fermi-LAT threshold. Thus, we 
may write 



dE^ 



signal dE. 



7 



prompt dE. 



7 



ICS 



and the total flux ^.y{Eth) as a function of the threshold energy Eth is 



Eth ^-^7 



%{Eth)= / dE,^^ . . (2) 



signal 



In the following, we discuss in detail our calculation of the prompt and the ICS components 
in the supersymmetric 7 signal. 

3.1.1 The prompt component 

The prompt component may be written as 



d%{E. 



1) 



dE.y 



1 1 ^dN. 



prompt AtT 2 ^ 



dN^ r r 

/ dQ / p^r{s,^|J))ds, (3) 



where we sum over the channels, X, of neutralino pair annihilation, both monochromatic and 

7 



continuum. The dN^/dE^ are the photon energy spectra produced per annihilation event. 
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for each channel with partial annihilation cross section (af)^, which we may separate as 
^ ^^^^^^ = '^^^^..^(E,-m^)+vaz,5\E, - - -^)) +Y. ^ ^) ^ ■ (4) 

The first and the second terms correspond to the monochromatic channels XX ~^ '^l/Zli 
respectively, for which we use the one-loop annihilation cross sections va^^ and vaz-f com- 
puted in [521 [53]. For the continuum component we calculate the dN^^^^/ dE^ using the 
partial cross sections (cf)^ for each of the two-body neutralino pair annihilation channels 
i = {qq ,r^T~ ,W~^W~ , ZZ , gg , Zh}, multiplied by the corresponding continuum spectra 
modelled using PYTHIA. 

We follow the standard procedure [5] for the halo factor in ([3[), defining a dimensionless 
average halo factor as an integral along light-of-sight (los) directions within a solid angle 

J{An)An = -^ [ dn [ p^{r{s,^))ds. (5) 
^&P& J An Jios 

where Rq is the solar distance from the GC and is local dark matter density. With these 
definitions, the prompt part in ([3]) takes the form 



dE^ 



prompt 2tx \2m 



\ ^ - ^ dN^ 



X ■ 



(6) 



X/ X 

We use the values Rq = 8.5 kpc and = 0.3 GeV/cm^. The integration along the light-of- 
sight direction can be performed using 

r'^{s,'ip) = Rl- 2s Rq cos ijj + s'^, (7) 

where ip is the angle with respect to the GC direction. Defining the function 

J{iIj) = —^ p\jRl-2sRQCOs^ + s^)ds, (8) 

J^qPq Jo ^ 

one can calculate the halo average in a cone within angle (f), that defines the size of the 
observation window around the GC, as 

J(Afi) An = 2n I J{ij)d{cos'ip) . (9) 

J cos 

This is the quantity introduced in ([5|), where we have used AQ = 27r(l — cos</)). The upper 
limit Smax in & is calculated in terms of the Milky Way halo size -Rmw 



= ^^/Rl^ - (i?osinV')2 + Rq cos^ (10) 
10 



although, as noted in [5l], the contribution of the integration beyond the scale radius 
~ 20 — 30 kpc is neghgible. 

The behavior of the galactic DM halo, especially in the inner region, is a very controversial 
topic. Although some N-body simulations indicate that the halo presents a highly cusped 
behavior towards its center [36l[55], other simulations [39l[56] as well as atomic hydrogen 
Hi observations on dwarf spiral galaxies [57] indicate shallower profiles. Recently, the Via 
Lactea 2 simulations [58] seem to support a cuspy profile like NFW, whereas the Aquarius 
DM simulation project [59] favors a different parameterization [38|l60]. that does not have 
such a cuspy character in the inner region, such as the Einasto profile. This unsettled 
issue of the cuspiness has significant impact on the theoretical predictions for the 7-ray flux 
emitted from the GC, and suggests a substantial theoretical uncertainty in calculating this. 
To address this uncertainty, we study three halo profiles that behave differently in the inner 
region of the Milky Way: NFW [36], Einasto [371138] and a simple isothermal profile [39] . 
The NFW behaves like at small distances, while the Einasto and isothermal profiles are 
both non-singular towards the galactic center. 

These profiles are defined by the density functions: 



We choose the following values of the model parameters: = 20 kpc, r'^ = 5 kpc and 
a = 0.17. The ps parameter, which has a different numerical value for each profile, is defined 
in such a way that Pq = p{Rq) = 0.3GeV/cm^. Assuming these numerical values, values 
of the halo factor J(Af2) AQ for three windows around the GC: 10, 7 and 5 degrees are 
tabulated in Table [1] The values in this Table are useful for understanding the relative sizes 
of the halo factor for the three halo profiles we study and the various windows. In the case of 
the the prompt component, the halo factor is simply multiplicative. In the case of the the ICS 
component, the halo factor is also, to a good approximation, multiplicative [H], but there 
are also other terms related to the photon background parameters and to the propagation 
of the lepton fluxes, etc., that are different for the various halo profiles and windows. 




Piso('^) 



(11) 



1 + (r/r^)2 ■ 
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Model 



10 deg 7 deg 5 deg 



NFW 
Einasto 
Isothermal 



10.51 7.90 5.95 
19.68 15.21 11.56 
1.21 0.62 0.32 



Table 1: Values of the halo factor J(Af2) AQ for various galactic dark matter halo profiles 
and values of the angle that defines the size of the observation window around the GC. 

3.1.2 The Inverse Compton Scattering component 

Neutralino pair annihilation produces energetic electron and positron (e^) fluxes mainly 
indirectly through the hadronization, fragmentation and decays of the primary annihilation 
products 0. The indirectly produced scatter on the ambient photon background, and 
through the ICS effect produce energetic photons. In order to calculate the ICS part we 
follow the tools and methods described in [TU|l7f]. 
The ICS flux is given by the equation 



where Eg is the injection energy and the halo function for the inverse Compton radiative 
process Iic{E:y, Es,ip) is defined as 



annihilation in the halo using PYTHIA, and the index i runs over the channels used for the 
photon prompt part. The functions V^q, the energy loss coefficient function b{E, r) and the 
generalized halo function I{E, Es,r{s,ip)) are given in |17] and references therein. 

The energy loss coefficient b{E, r) depends on the profile of the magnetic field in the 
galactic plane, as described in |17l|6l]. Recent estimates suggest that the magnitude of the 
magnetic field towards the GC can be larger by a factor of ten or more than in the solar 

■^The fluxes produced directly through the channels 2% — e+e^ are negligible, because this process is 
P-wave suppressed and proportional to the small electron mass. 




(12) 




(13) 



We calculate the primary electron/positron fluxes dN^^/dE produced by neutralino pair 
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neighbourhood [62]. We have verified that a magnetic field that is stronger by a factor of 10 
would result in a reduction of the ICS fiux almost by a factor of two. This would affect the 
total 7-ray flux at low photon energy, where the ICS is dominant. However, since the most 
important part of the CMSSM 7 signal is at higher energies, its effect on the global that 
we calculate in the next Section is small, of the order of 10%. A similar remark applies to 
the possibility of energy losses due to interactions with gas in the GC [63] . 

The sum over the index a includes three photon backgrounds: the cosmic microwave 
background (CMB), the starlight in the galactic plane and the infrared radiation due to the 
rescattering of the starlight by dust. To calculate the ICS part, we follow the semi- analytic 
method described in [U], see also [TOlIM] . which yields results similar to numerical methods 
such as GALPROP [65]. In particular, for the propagation parameters there are three 
commonly used models, called MIN, MED and MAX [61], that correspond to the minimal, 
median and maximal primary positron fluxes compatible with data on the boron to carbon 
ratio, B/C [66] • In this paper we choose the MED set of parameters, but we have checked 
that our flndings do not depend signiflcantly on this choice. Specifically, scanning through 
the various parameter sets, we find that the value of the ICS fiux can be changed by at most 
15%. In the low-energy range, this uncertainty affects analogously the total 7-ray signal, 
but its effect in the statistical analysis in Section H] is unnoticeable. 

3.2 Examples within the CMSSM 

Using the above treatments of the prompt and the ICS parts, we can evaluate the total flux 
in (III) and ([2]) at any point of the parameter space of the CMSSM. In Fig. HJ we display 
examples of the differential flux E^d^^/dE^ as a function of the photon energy E^ in a 
7-degree window around the GC. The blue curves are the prompt parts and the green curves 
correspond to the ICS contributions. We use the NEW proflle in each case, and the four 
panels in Eig. H] correspond to the four benchmark points C,E,L and M introduced in [31] 
and updated in [67]. These points belong to different characteristic regions of the CMSSM 
model that are cosmologically acceptable. 

Referring to the left panel of Eig[T]for tan /3 = 10, point C represents the low-mass region 
of the coannihilation strip for tan /3 = 10 with (mi/2,''Tio) = (400,96) GeV, and point E is 
in the focus-point region with {mi/2,mQ) = (300,2003) GeV (outside of the range shown 
in Eig [T^). Point L is in the coannihilation region at (mi/2,mo) = (450,312) GeV and 
tan P = 50, and referring to the right panel of Eig [1] for large tan /3, point M is in the rapid- 
annihilation funnel region at large [7711/2,1^0 = (1075,1045) GeV with tan/3 = 55. The LSP 
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Figure 4: The differential flux d^^/dE^ as a function of the photon energy E^ for a 7 deg 
window around the GC. The blue curve is the prompt part and the green curve is the ICS 
flux. We use here the NFW profile, and the four panels correspond to the four benchmark 
CMSSM points C,E,L and M, which are described in the text. 



masses at the points C,E,L and M are = 165, 117, 193 and 474 GeV, respectively. Since 
point C is in the coannihilation region, it has a small neutralino annihilation cross section, 
as we saw in Fig. [2], and hence the photon flux for point C is weaker than for the other 
benchmark points. The focus-point region point E has a stronger flux, because it has a 
larger annihilation cross section and small m^. The points L and M have comparable fluxes, 
but that of point M is relatively smaller, because there is larger ^. 

In Fig. [5l we plot the fluxes ^^{Eth) in a 7-degree window around the GC, for the same 
benchmark points as functions of the threshold energy Eth., as defined in ([2]). Once again, 
the green curves represent the ICS parts, the cyan curves represent the monochromatic parts 



""It is useful to note here that both the prompt and ICS fluxes are, in general, inversely proportional to 
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E,h [GeV] E,^, [GeV] 

Figure 5: The total flux ^^{Eth) as a function of the threshold energy Eth for a 7 deg window 
around the GC. The green curves represent the ICS part, the cyan curves represent the 
monochromatic part, and the blue curves the continuum part. We use here the NFW profile, 
and the four panels correspond to the four benchmark CMSSM points C,E,L and M, as in 
Ftg.^ 

and the blue curves the continuum parts. We see in these figures that the monochromatic 
component yields only a very small contribution in every case, as was discussed in the previ- 
ous Section, and is also inversely proportional to m^. Another remark concerns the relative 
magnitudes of the ICS and prompt components. Since the ICS component results from 
background photons energized by interaction with energetic charged leptons, it dominates 
for photon energies much smaller than m^. In particular, for LSP masses up to a few hundred 
GeV as along the benchmark strips we study, the ICS contribution becomes subdominant 
above the range -E^ > 1 — 3 GeV, as can be seen in Fig. HI 
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3.3 The background and Fermi-LAT data 

The background has three main components: diffuse galactic emission (DGE), an isotropic 
extragalactic contribution and resolved point sources. DGE is produced essentially by the 
interaction of cosmic-ray nucleons and electrons with the interstellar medium. Nucleon inter- 
actions produce photons via vr^ decays, whereas electrons produce DGE via bremsstrahlung 
and ICS. For the determination of the DGE background component, a conventional model 
was developed [M] using the hypothesis that the cosmic-ray spectra in the Galaxy can be 
normalized to solar-system measurements. This model did not reproduce well the EGRET 
data, especially for ~ C(GeV), resulting in the apparent EGRET excess [69]. An up- 
dated version of the model that was able to describe the data was developed in [70] , and the 
recent Fermi-LAT measurements on the diffuse background component [12] are well repro- 
duced by the updated version of the conventional model. This is the diffuse model that we 
use in our analysis. As was pointed out in [70], to estimate the errors in this model is a 
complicated task, but based on the quality of its fit to the data Fermi-LAT, we can say that 
is about 20% - 25% for photon energies up to 100 GeV. 

The isotropic diffuse emission or isotropic gamma-ray background (IGRB) is much fainter 
and, although the term extragalactic gamma-ray background is usually applied, its extra- 
galactic origin is not universally accepted. Various astrophysical objects can contribute to 
this emission, such as active galactic nuclei, galaxy clusters, blazars, star-forming galax- 
ies, ultra-high-energy cosmic rays, etc. [7T]- The Fermi-LAT collaboration measurement of 
the IGRB [15] is consistent with a power law, somewhat softer than the previous EGRET 
analogous measurement [72]. In any case, the isotropic component is subdominant in our 
analysis. 

On the other hand, resolved point sources (RPS) provide an important part of the pho- 
ton background from the direction of the GC. We use the first 11-month Fermi-LAT cata- 
logue [H], which contains 1451 point sources modelled in the energy range 100 MeV to 100 
GeV. The gamma flux of each point source is taken to obey a simple power law. 

In order to evaluate the possible constraints imposed on the CMSSM parameter space by 
the latest Fermi-LAT data [13], we estimate the background using the Fermi-LAT Science 
Tools [73]. In our analysis, we use data collected by Fermi-LAT between Aug 4, 2008 and 
April 29, 2011, making a selection based on the recommendations of the collaboration. We 
focus our analysis in the energy range 300 MeV to 300 GeV, dividing it in 25 bins spaced 
evenly on a logarithmic scale. We have studied various windows around the GC in the 
range 1-10 degrees, but we choose as the basic region-of-interest (ROI) of our analysis the 
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Figure 6: The various background components in counts per MeV for a 7 deg window around 
the GC in the range 300 MeV to 300 GeV, estimated using the Fermi-L AT Science Tools J75y . 
We also display the Fermi-L AT data ITB^ collected during the period described in the text. 



7-degree window centered at the position of the brightest source in the GC: RA = 266.46", 
Dec = —28.97°, as in [16]. Note however, that our conclusions are not overly sensitive to 
this choice of window size. 

Fig. |6] displays the various background components obtained using the Science Tools [73] . 
along with the data collected during the aforementioned period. The error bars attached 
to the data represent the purely statistical errors. We have performed a binned likelihood 
analysis using the gtlike tool, including in the background model file used in the fit 1) the 
galactic diffuse model, 2) the isotropic spectral template, and 3) the point sources, all as 
provided by the collaboration. Concerning the point sources, we have included additional 
RPS from the vicinity of the ROI, that possibly can contribute to the observed counts. 

In Fig. [7] we plot the residuals in each of the 25 bins in the range 300 MeV to 300 
GeV, i.e., the difference (counts — model) /model, where in this case "model" includes just 
the background sources. As discussed earlier, the background models used to describe the 
collected data are able to reproduce them with an accuracy on the order of 25%, especially 
for photon energies up to 100 GeV. This was expected, since the background models are not 
adequate for energies larger that this |T2l[T3]. The shaded area represents the uncertainty 
in the effective area of the detector. As discussed in [16], the systematic uncertainty in the 
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Figure 7: The residuals in each bin of the available Fermi-LAT data JT^ used in Fig. 0, 
for a 7 deg window around the GC in the range 300 MeV to 300 GeV, as compared to the 
estimated background model. The shaded area represents the uncertainty in the effective area 
of the detector as discussed in the text. 

effective area of the LAT is currently estimated as 10% at 100 MeV, decreasing to 5% at 560 
MeV, and then increasing to 20% at 10 GeV. As suggested by in [12], we assume that this 
uncertainty propagates to the model predictions. With the systematic uncertainty included, 
we find = 31.1 for 25 degrees of freedom, corresponding to a p- value of 19%. 

4 Experimental Perspectives 

Using the results of the previous Section, we now compare the signals expected in CMSSM 
scenarios with the background, under various hypotheses, again using 25 logarithmic bins in 
the range 300 MeV to 300 GeV. We define the function as 



where di is the number of the data counts per bin, bi is the expected background as calculated 
by Science Tools [73], and Si is the signal due to the dark matter halo for the corresponding 
bin. For the uncertainty, we have included both the statistical uncertainty in the data count 
as well as the systematic uncertainty in the effective area, (Tea- These have been added in 
quadrature, so that erf = di + a^^. In Fig. |8l we present the current Fermi-LAT data set [13] 
as blue dots with statistical error bars, the background estimate (red curve) and with the 
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Figure 8: The current Fermi-LAT data set fT^ . in counts per bin, as a function of photon 
energy in MeV, including the statistical error (blue points) for a 7-degree window around 
the GC. The red line is the corresponding estimated background. We also display the signals 
expected for the three halo profiles we are using: NFW, Einasto and isothermal and the four 
panels correspond to the four benchmark CMSSM points C,E,L and M, as in Fig. [7[ 

signal corresponding to the three DM halo profiles discussed previously: NFW (green), 
Einasto (purple) and simple isothermal (cyan). The four panels of these figures correspond 
to the four CMSSM benchmark points C, E, L, M discussed in the previous Section, so the 
relative sizes of the signal can easily be understood. The coannihilation point C produces 
the weakest signal, as expected, and the focus-point benchmark E gives the strongest signal. 
Points L and M yield signals that are almost comparable, again as expected from Fig. |H 
Had we chosen a point from the right side of the funnel, the signal would have been about 
a factor of 2 higher than shown for point M. 

We start with a nominal implementation of the current experimental situation, that uses 
the present Fermi-LAT data set [13], assumes the current understanding of the background 
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and systematic errors, and uses an NFW density profile. As a preliminary to examining the 
signals along the WMAP strips shown in Fig. [H we first consider a more global view across 
an (mi/2,mo) plane in the CMSSM. As described in the previous section, in the absence of 
any supersymmetric contribution, our modeling of the astrophysical background give a value 
of = 31.1 for 25 degrees of freedom. In Fig. [9l we show for illustration (in yellow) contours 
of constant Ax^ across the CMSSM {mi/2,mo) plane for tan/3 = 55 exhibiting contours 
of Ax^ = 0,1 and 4 relative to the value of without supersymmetry Q Note that there is 
a small region at low mi/2, in the region already excluded by 6 — )■ S7, where there is a small 
but insignificant improvement to x^. As one can see in the left panel of Fig. |9l where an 
NFW profile has been assumed, a modest constraint is imposed by the present Fermi data. 
In addition to the area at small mi/2 that is already excluded by supersymmetry searches at 
LEP, the Ax^ = 1 contour runs between the branches of the rapid-annihilation funnel, and 
disfavours the right branch at the 68% CL. The Ax^ = 4 contour runs deep inside the funnel 
where the relic density is below the WMAP range. In the right panel, where the Einasto 
profile has been assumed, the disfavoured region is somewhat larger, as expected from the 
values of the halo factors in the Table [H On the other hand, based on these values, the 
opposite is expected for the for the isothermal model. 

We now focus on the WMAP strips for tan/3 = 10 and 55. The results of our nominal 
implementation of the current and possible future Fermi constraints are shown in Fig. [10] 
for the coannihilation and funnel strips (upper panels) and tan/3 = 10 (left panels) and 
tan/3 = 55 (right panels). As a function of mi/2 along each of these lines, we display the 
X^ function for the background alone, which is a constant 31.1 for 25 degrees of freedom, 
and the corresponding x^ function for the combination of the background with the signal 
calculated as a function of mi/2 in each scenario, for the current Fermi-LAT data [I3]. Also 
shown are the possible results from future 5- and 10- year data sets. The baseline x^ is 
correspondingly higher due to the expected improvement in the statistical uncertainty in the 
data, assuming no improvement in the systematic uncertainty associated with the effective 
area. The horizontal lines assuming no supersymmetry are not shown for these hypothetical 



^We do not show analogous results for tan (3 — 10, since the results arc unpromising, as we see below 
when we restrict our attention to the WMAP strips. 

^We recall that the relic LSP density is too large in the interior of the plane between the WMAP strips, 
and too small in the outer regions between the WMAP strips and the charged-LSP and electroweak symmetry 
breaking boundaries. The interior region is forbidden in the i?-conserving scenario discussed here, whereas 
the outer regions could be acceptable in the presence of some other component of dark matter. However, 
in this case the annihilation rate should be reduced by the square of the relic LSP density relative to the 
WMAP dark matter density. Only the intersections of the Ax^ contours with the WMAP strips should be 
taken literally. 
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Figure 9: Contours of Ax^ = 0, 1,4 in the plane /or tan /3 = 55 (shown in yellow). 

In the left panel, we assume an NFW profile for dark matter in the core of the galaxy, whereas 
in the right panel we assume an Einasto profile. Conventions for the other curves and shaded 
regions are as for Fig. [11 



data sets. In each case, the results assuming an NFW profile are shown as solid lines, and 
results assuming an Einasto profile are shown as dashed lines. 

We see immediately that the calculated signal with the current Fermi-LAT data set [13] 
has very little effect on the overall function along the coannihilation strip for tan /3 = 10 
(upper left panel), whereas somewhat larger effects are visible along the WMAP strips shown 
in the other panels. However, in all cases except along the second funnel strip for tan /3 = 55 
the change in the function is < 1 in the case of the NFW profile (solid blue lines), 
and hence not significant. In general, adding the supersymmetric signal worsens slightly 
the quality of the fit, though there is a slight improvement along the focus-point strip for 
tan /3 = 10 for 200 GeV < mi/2 < 400 GeV. The (red) dashed curves in Fig. [10] show 
the resulting functions assuming an Einasto dark matter profile. While this does yield 
somewhat larger effects, we see that Ax^ < 3 even along the right branch of the funnel for 
tan /3 = 55, and a 95% CL exclusion is not possible. 

Hence, the most that one could hope for along most of the WMAP strips studied is some 
inconclusive level of disfavour, except in the limited region of small mi/2 for tan /3 = 10 where 
some improvement in the overall x^ of a global fit might eventually be possible. Significant 
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Figure 10: The functions along the CMSSM WMAP strips as functions of mi/2 for 
tan /3 = 10 (left panels) and tan/3 = 55 (right panels), in the coannihilation and funnel 
regions (upper panels) and in the focus-point region (lower panels). In each panel, we display 
the function for the background alone, a horizontal line at = 31.1, and the x"^ function 
obtained by adding the calculated x " X annihilation signal in the current (approximately 2 
1/2 year) Fermi data sample and in projected 5- and 10-year data sets. Solid (blue) curves 
are based on an NFW profile, while dashed (red) curves are based on an Einasto profile. 
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Figure 11: The functions as functions of 1711/2 along the CMS SM focus-point WMAP 
strip for tan /9 = 55 assuming a factor of 2, 3 or 4 reduction in the systematic error (left 
panels) or a negligible systematic error (right panels) , assuming in all cases that an improved 
estimate of the background brings it to ±lcr from the data. The upper (lower) panels are for 
the coannihilation/ funnel strip (focus-point strip). 
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Figure 12: As in Fig. [771 but for tan/3 = 10. 



improvements in our current understandings of the background and/or systematic errors 
would be required before the scenarios studied could be constrained significantly. 

Examples of the effects of possible improvements in understanding the background and 
systematic errors are shown in Fig. [HI in the cases of the coannihilation/funnel strip (upper 
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panels) and focus-point strip (lower panels) WMAP strips with tan /3 = 55, assuming an 
NFW profile. The functions in the absence of any supersymmetric signal are shown as 
horizontal black lines, and the effects of including a signal are shown as blue lines. In all of 
these panels, we have assumed that improvements to our understanding of the background 
bring it to lie within 1 a of the data, with the sign of the discrepancy chosen randomly in 
each bin. In the left panels, we also assume reductions in the current systematic error in the 
effective area by a factors of 2, 3 and 4 in 0"^^, whereas in the right panels the systematic 
error is assumed to be negligible B We see immediately in the left panels that reducing 
the systematic error by a factor of 2 does not improve greatly the discrimination between 
the background-only and signal -|- background hypotheses. Even a reduction by a factor 
of 4 does not lead to a 95% CL exclusion limit, whereas we see in the right panels that 
a reduction in the systematic error to a negligible level might provide such discrimination. 
Depending on the evolution of the understanding of the background, adding a CMSSM 
might affect significantly the fit for small mi/2- With the present background estimate, 5 
or 10 years of Fermi data might provide evidence for the CMSSM at small mi/2, albeit in 
the region disfavoured by 6 — )■ 57, or possibly evidence against the CMSSM in the region 
mi/2 > 500 GeV where accelerator constraints are currently irrelevant. The corresponding 
figures for tan /3 = 10 with improved systematics are shown in Fig. [121 In the case of the 
coannihilation strip for tan /3 = 10, even the absence of systematic errors would not allow 
significant information to be obtained, though some information might be obtainable about 
the focus-point strip. For instance, using the Einasto profile instead of NFW, the of the 
last panel in the Fig. [12] is about 3 times bigger, enabling more definite conclusions to be 
drawn. 



5 Summary and Prospects 

We have studied in this paper the potential of searches by the Fermi-LAT detector [13] 
for 7 rays from the Galactic centre produced by the annihilations of neutralino LSPs in 
the CMSSM. We have shown predictions for benchmark points, {'mi/2,Tno) planes, and the 
concentrated on the WMAP strips for tan/3 = 10 and 55. We have found that the 7-ray 
signal would be very difficult to detect along the coannihilation strip for tan/3 = 10, as 
expected from the relatively small XX annihilation cross section along this strip. The focus- 

^For large even a small signal improves slightly the fit, in comparison to the background case only. So 
we see the of the signal is smaller than the of the background for large mi/2- This improvement, for 
the case of the reduced systematic error, is due to the contribution of the last bins that contain energetic 
photons with > 150 GcV. 
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point strip for tan/3 = 10 and both the coannihilation/funnel and focus-point strips for 
tan /3 = 55 have larger annihilation cross sections and offer better prospects, though even in 
these cases the present data do not offer good prospects for discriminating between models, 
because of uncertainties in the background estimates and the current systematic errors of 
the Fermi-LAT detector. 

Future Fermi-LAT data sets might offer some information about parameters of the 
CMSSM, also in regions of parameter space not accessible to accelerator experiments. How- 
ever, this would require several more years of data and a substantial reduction in the present 
systematic error. This is therefore a priority for optimizing the prospects for future Fermi- 
LAT constraints on the CMSSM from searches for 7 rays from the Galactic centre. Quan- 
titative understanding of the constraints imposed by Fermi-LAT results would also require 
advances in our understanding of the dark matter density profile in the core of the Galaxy, 
since we see considerable differences between results assuming NFW and Einasto profiles. 

The Fermi-LAT detector also has data on 7 rays from many other astrophysical sources, 
including the Galactic bulge, dwarf galaxies and diffuse emissions. We therefore plan to 
extend our analysis to these cases. However, it is clear that it will always be more difficult to 
pin down the CMSSM than some other supersymmetric models via searches for energetic 7 
rays from astrophysical sources, and the same remark applies to other astrophysical strategies 
to look for neutralino annihilations. 
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